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Abstract 

We present the first quenched lattice QCD study of the form factors relevant for 
the hyperon semileptonic decay E~ n (. v. The momentum dependence of both 
axial and vector form factors is investigated and the values of all the form factors 
at zero-momentum transfer are presented. Following the same strategy already ap- 
plied to the decay tt~£ v, the SU(3)-breaking corrections to the vector form 
factor at zero-momentum transfer, /i(0), are determined with great statistical accu- 
racy in the regime of the simulated quark masses, which correspond to pion masses 
above k, 0.7 GeV. Besides /i(0) also the axial to vector ratio 5i(0)//i(0), which 
is relevant for the extraction of the CKM matrix element Vus^ is determined with 
significant accuracy. Due to the heavy masses involved, a polynomial extrapolation, 
which does not include the effects of meson loops, is performed down to the physical 
quark masses, obtaining /i(0) = -0.948 ± 0.029 and gi(0)//i(0) = -0.287 ± 0.052, 
where the uncertainties do not include the quenching effect. Adding a recent next- 
to-leading order determination of chiral loops, calculated within the Heavy Baryon 
Chiral Perturbation Theory in the approximation of neglecting the decuplet contri- 
bution, we obtain /i(0) = — 0.988±0.029iattico±0.040HBChPT- Our findings indicate 
that SU(3)-breaking corrections are moderate on both /i(0) and .gi(0). They also 
favor the experimental scenario in which the weak electricity form factor, (?2(0), is 
large and positive, and correspondingly the value of |<?i(0)//i(0)| is reduced with 
respect to the one obtained with the conventional assumption g2{<f') — based on 
exact SU(3) symmetry. 



1 Introduction 



Recently it has been shown that SU(3)-break;ing corrections to the K ^ it vector form 
factor can be determined from lattice simulations with great precision [1]. The approach 
has allowed to reach the percent level of accuracy in the extraction of Vus from Kf^^, 
decays, thus stimulating new unquenched lattice studies to reduce the systematic un- 
certainty [2]. An independent way to extract Vus is provided by hyperon semileptonic 
decays, and Ref. [3] has shown that it is possible to extract the product \Vus ■ /i(0)| at 
the percent level from hyperon experiments, where /i(0) is the vector form factor (f.f.) 
at zero-momentum transfer. The Ademollo-Gatto (AG) theorem [4] protects /i(0) from 
first-order SU(3)-breaking corrections that are thus suppressed. Experiments seem to 
be consistent with negligible channel-dependent SU(3) corrections in this f.f. [3], but 
they cannot exclude sizable (i.e. larger than percent), channel- independent effects in the 
extraction of Vus- Model dependent estimates based on quark models, and chiral 

expansions give different results (see e.g. Ref. [5]), so that a lattice QCD determination 
of /i(0), as well as of all the other vector and axial f.f.'s which are not AG protected, is 
of great interest. 

The aim of the present work is to investigate SU(3)-breaking corrections to both 
vector and axial form factors relevant for the — n transition, completing and gen- 
eralizing in this way the preliminary lattice study of Ref. [6]. Though our simulation 
is carried out in the quenched approximation, our results represent the first attempt to 
evaluate hyperon f.f.'s using a non-perturbative method based only on QCD. 

We first show that it is possible to determine the SU(3)-breaking corrections to /i(0) 
on the lattice following the method of Ref. [1], which is based on the following three main 
steps: i) evaluation of the scalar form factor /o('Z^) at = Q'jnax = (-^s — -^n)^; h) study 
of the momentum dependence of fo{q^) to extrapolate /o(Q'max) down to /o(0) = /i(0); 
and iii) extrapolation in the quark masses down to the physical point. The latter step 
is one of the main sources of uncertainty in the quenched lattice calculations of the 
vector f.f., since the AG theorem makes this quantity dominated by meson loops. We 
make use of the recent, next-to-leading order (NLO) calculation of the chiral corrections 
to /i(0) performed in Ref. [7] within the Heavy Baryon Chiral Perturbation Theory 
(HBChPT). As in the meson sector, the AG theorem prevents the contribution from 
local counter-terms up to O(p^) and makes the NLO corrections finite and free from 
unknown low-energy parameters. However, the convergence of the chiral expansion turns 
out to be rather poor and the inclusion of the decuplet contribution appears to spoil the 
expansion itself [7]. Thus, the third step of the procedure of Ref. [1], i.e. the correction 
of the leading quenched chiral logs with the full QCD ones, is not possible in the hyperon 
case. In this respect we note that the use of HBChPT in conjunction with lattice QCD is 
likely to represent a serious limitation for achieving a precise determination of /i (0) also 
in future (unquenched) lattice calculations, until it will be possible to simulate on the 
lattice very light quark masses, close to the physical ones, and sufficiently large volumes 
[8]. 

Due to the rather heavy masses involved in our simulation, it is unlikely that meson 
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loops can contribute significantly at the simulated quark masses. Thus a polynomial 
extrapolation, which does not include the effects of meson loops, is performed down to 
the physical quark masses, leading to /i(0) = —0.948 ±0.029. This result for the SU(3)- 
breaking corrections, corresponding to /i(0) + 1 = (5.2 ± 2.9)%, represents mainly an 
estimate of the local terms of the chiral expansion and appears to be opposite in both 
sign and size to the NLO estimate of chiral loops, (—4 ± 4)%, calculated in Ref. [7] in 
the approximation of neglecting the decuplet contribution and assuming a 100% overall 
uncertainty. The two contributions should be added, leading to our final estimate: 
/i(0) = —0.988 ± 0.029iattice i 0-040HBChPT) where the uncertainties do not include the 
quenching effect. 

Second, we investigate the momentum dependence of the other two vector f.f.'s, the 
weak magnetism f2{q^) and the induced scalar fsiq^), as well as of the three axial f.f.'s, 
the axial- vector gi{q'^), the weak electricity 52 (q^) and the induced pseudoscalar (73(5^)^. 

The ratio gi{q^) / fi{q^) can be determined with quite good statistical accuracy di- 
rectly from the ratio of appropriate axial and vector three-point correlation functions. 
We show that such a ratio can be evaluated on the lattice at q^ = q^g^ (i.e. with hyper- 
ons at rest) within ~ 10% accuracy and that the extrapolation down to = docs not 
modify significantly that level of precision. Adopting a polynomial extrapolation in the 
quark masses we get gi{0) / fi{0) = —0.287 it 0.052, which is consistent with the value 
ffi(0)//i(0) = -0.340 ± 0.017 adopted in the recent analysis of Ref. [3]. The study of 
the degenerate transitions also allows to determine the value of (7i(0)//i(0) directly in 
the SU(3) limit; we get [gi{0)/ fi{0)]su{3) = -0.269 ± 0.047. Compared with the above 
estimate, it implies that the SU(3)-breaking corrections on yi(0)//i(0) are moderate, 
though this ratio is not protected by the AG theorem. Our finding is in qualitative 
agreement with the exact SU(3)-symmetry assumption of the Cabibbo model [9]. 

Contrary to fi{q^) and gi{q'^), the other f.f.'s cannot be evaluated at = g^ax 
and consequently the extrapolation of the lattice data to q^ = is affected by larger 
uncertainties. In the cases of the weak magnetism and of the induced pseudoscalar 
f.f.'s, whose matrix elements do not vanish in the SU(3) limit, we obtain /2(0)//i(0) = 
— 1.52 lb 0.81 and 53(0)7/1(0) = 6.1 it 3.3 at the physical point. The central values 
agree well with the experimental result /2(0)//i(0) = —1.71 ± 0.12stat. ± 0.23syst. from 
Ref. [10], as well as with the SU(3)-breaking analysis of Ref. [11], and with the value 
53(0)7/1(0) = 5.5 ± 0.9, obtained using the axial Ward Identity and the generalized 
Goldberger-Treiman relation [12], which relates 53(0) with 51 (0) in the chiral limit. 

For the weak electricity 52 ((Z^) and the induced scalar fsiq^) f.f.'s a non- vanishing re- 
sult, which is entirely due to SU(3)-breaking corrections, is found, namely 52(0)7/i(0) = 
0.63±0.26 and /3(0)7/i(0) = -0.42±0.22. Note that experiments carried out with polar- 
ized S~ hyperons [10] have determined the f.f. combination |5(i(0)7/i(0) — 0.133 52(0)7/1(0)1 
to be equal to 0.327 ± 0.007stat. =t O.OlOgyst.- Our result is 0.37 ± 0.08, which means that 
the lattice results for both 5i(0)//i(0) and 52(0)7/i(0) are nicely consistent with the 
experimental data on the S~ — > n transition. Our findings favor the scenario in which 

^Our definitions of the f.f.'s /2, f,i, (?2 and gs [sec Eq. (1)] differ by the ones adopted in Ref. [3] by a 
factor (Me -|- M„)/Me, equal to ~ 1.785 at the physical hyperon masses. 
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52(0)//i(0) is large and positive, and correspondingly the value of |(7i(0)//i(0)| is re- 
duced with respect to the one obtained with the conventional assumption g2{(f') = 
(done in Ref. [3]) based on exact SU(3) symmetry. Such a scenario is slightly preferred 
also by experimental data (see Ref. [10]). 

The plan of the paper is as follows. In Section 2 we introduce the notation and give 
some details about the lattice simulation. Section 3 is devoted to the extraction of fo{q^) 
at = (/max) to its extrapolation down io q'^ = and to the polynomial extrapolation in 
the quark masses down to their physical values, without taking into account the effects of 
meson loops. The estimate of these effects, based on HBChPT at the NLO, is described 
in Section 4, where the problematic related to the convergence of the chiral expansion 
and to the decuplet contribution is briefly illustrated. In Section 5 we present our results 
for the ratio 5i(0)//i(0), while those for the other vector and axial f.f.'s are collected in 
Section 6. Finally our conclusions are given in Section 7. 



2 Notation and Lattice Details 

In this study we restrict our attention to the n I v decay. We are interested in 

the hadronic matrix clement of the weak {V — A) current, M = {n\u'y^{l — 75),s|E~), 
which can be conveniently expressed in Minkowski space in terms of f.f.'s and external 
spinors as 



M = u„(p') 7''/i(9')-^ 



+ 



52 (g') + 



53(g') 



75 \ UT,{p) 



= u„(p'){0(^ + 0^}^s(p) 



(1) 



with q = p — p'- A detailed discussion of the properties of the six f.f.'s in Eq. (1) can 
be found, e.g., in Ref. [3]. We also introduce the scalar f.f. fo{q'^) from the divergence of 
the vector weak current {n\dn Vf^\T,~) = (Ms + M„)/o(g^). It is related to /i(g^) and 
f3{q') by 



/o(g')=/i(g') + 



(2) 



n 



and at gr^ = it coincides with the quantity of interest /i(0). Note that for the S" 
transition /i(0) is normalized as /i(0) = —1 in the SU(3) limit. 

In order to access the matrix element of Eq. (1) on the lattice, one can consider the 
following two- and three-point correlation functions in euclidean space 



G^^-\t,p} 



^(J^^/")(^,x)J^(")(0,0))e-^^-^ 



7'7 



-i(p-p')-x -ip'-y 



[V^^it,,ty,p,p')]^,^ = {J:;,{ty,y)Vi,it,,x)J^{0,0))e 
[A^''{tx,ty,p,p')L,^ = (j;,(t„y)A^(t.,x)J^(0,0))e-(^^-^^')- e-'P'-y , 



(3) 

(4) 
(5) 
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where Jy ' are local interpolating operators for the neutron and the hyperon, which 
we choose to be 

with latin (greek) symbols referring to color (Dirac) indices. The operators in Eq. (6) 
are related to the external spinors us,n (cf- Eq. (1)) via the relation 



(0|jf")(0,0)|S(n)(p,a)) = [uj,^n)ip, '^)]^ , (7) 

where a refers to the polarization of the S(n) baryon. 

In Eqs. (4)-(5) and are the renormalized, 0(a)-improved lattice weak vector 
and axial-vector current respectively: 

Vf" = Zv(l + by^^I^l±^] {u^^^s + acvd, ua'^'^s) , (8) 



A^^ = Za{1 + 6^^!!!£±^ ) ^ur^^s + acAd^ u^^s) , (9) 



where Zy and Za are the vector and axial- vector renormalization constants, by, bA, cy 
and CA are 0(a)-improvement coefficients [13] and the subscript I refers to the light u 
or d quarks, which we consider to be degenerate in mass. 

Taking the large-time limits tx, (ty — tx) — >■ oo and using (7) one can rewrite Eqs. (3)- 
(4) as follows 



(10) 



'^^"A (oivf^^l , (11) 



'^'^""^^ iOiuC^) , (12) 



where = \J Ml + Ey, = y M| + jpj^ and O^Y {O^) is the euclidean version of 

the vector (axial-vector) contribution to Eq. (1). From Eqs. (10) and (11) it follows 

A^^{tx,ty,p,p')y^ 1 [(zi^' + M„)0;^(ij^ + Ms)]y7 ^^^^ 



G^{tx,p)^'JG^{ty-tx,p')y^*-^*y-*^^--'^/Z;;Z^ {i]^' + Mn)Yj (ij^ + Ms) 
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Let us consider the following kinematics 



Ml + \<f(^,q)^{E,,q), p' = {MrJ) (15) 

and study the matrix elements of the weak vector current. In order to minimize the 
number of calculated three-point correlation functions we consider only one pair of values 
of the Dirac indices, namely 7 = 7' = 0. Then we define the following quantities 



2Eq^/z;^ {y^'i^x^h^ 0)00) 

E, + Ms G^{t,, q)oo G-{ty - t,, 0)00 
2Eg^/Z;;Z^ I™ {ikV^^'itx, ty, q, 0)oo) 



Wi^-t t) ^ Re {Vi'''itx,ty, q, 0)00) 

' " ' qi G^{t,,q)oo G-{ty-t,,6U ' 

which in terms of the three vector form factors /i,2,3 read as 

Wstf) = /,(«") +/2(«') , (17) 

where Wi{q^) = lim^^^ _ ^ 00 Wi{q^;tx,ty) (i = 1, 2, 3). Inverting the above equations 
one gets 

/.(,^) ^ ^{W^i^) - ^TT.tf ) - , ^J^;^^ , ^,) W3fa-)}, ,18) 
/.(,^) ^ ^{-W.(,V|^W.tf) + ^W3b^)}. (19) 

with M = (Ms + Mn)/2Mn. 

Equations (16)-(20) provide the standard procedure for measuring f.f.'s on the lattice. 
Typically the accuracy that can be reached in this way is not better than 10-^20%, which 
is clearly not sufficient for measuring SU(3)-breaking effects in /i(0) at the percent level. 
In the next Section we describe the procedure to get both the scalar form factor (2) at 

= g^ax ^'^d the ratios f2{q^)/fi{q^) and h{<f)/ fi{q^) with quite small statistical 
fluctuations. 
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In the case of the axial- vector weak current, choosing always 7 = 7' = 0, we define 
the following quantities 



lE^^fZ^ ty. g, 0)00) 

+ Ms g-)oo G-{ty - t,, 0), 



' 9193 G^(tx,g1ooG-(t^-t,,0)oo 

where i^'^ = A^"^ + {Af^-fqi + Af"/g2) • (1^1^ - ql) /^Qs- In terms of the three axial- 
vector form factors 51,2,3 one has 

777(^)/ 2\ / 2\ Eq-Mn 2\ , Eg - 2\ 

= ^'^"^ ^ - M^TWJ'^' ^ + M^TWJ'^'^ ^ ^ 

w(^)/' 2\ f 2\ Eq + Mj^ 2\ , Eq-Mn 2^ 

W^i^^(9') = 92{q') - 93{q') , (22) 

where I^1'^^(q'^) = limj^, (t^ - 1,) ^ 00 t^;, t^^) (z = 1,2,3). Inverting the above 

equations one gets 

[(M„-Ms)(£;g-M„)-(M„ + Ms)(£^g-Ms)]}, (23) 
92iq^) = ^f{w^^V)-Wi'\q^) + ^;^Wi%% (24) 

gsiq') = M{w[''\q')-Wi''\q')-Wi''\q')} . (25) 

Wc have generated 240 quenched gauge field configurations on a 24^ X 56 lattice 
at (3 = 6.20 (corresponding to an inverse lattice spacing equal to ar^ ~ 2.6 GcV), 
with the plaquette gauge action. We have used the non-perturbatively 0(a)-improved 
Wilson fermions with csw = 1-614 [14] and chosen quark masses corresponding to four 
values of the hopping parameters, namely k G {0.1336,0.1340,0.1343,0.1345}. Using 
S and n baryons with quark content [kgkiki) and {k^kiki) respectively, twelve different 
E -H' n vector (axial) correlation functions V^" (^^") have been computed, using both 
kg < ke and kg > ke, corresponding to the cases in which the S (neutron) is heavier 
than the neutron(S). In addition, using the same combinations of quark masses, also 
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the three-point n ^ E correlation functions V^''^ have been calculated. Finally, 

twelve clastic, non-dcgcncratc V^'^ (^^^) and four elastic, fully degenerate V^^ 
three-point functions have been evaluated. 

The simulated quark masses are approximately in the range 1 ^ 1.5 x m^, where 
TTT-s is the strange quark mass, and correspond to pseudoscalar meson masses in the 
interval ~ 0.70 1 GeV and to baryon masses in the range ~ 1.5 1.8 GeV. Though 
the simulated meson masses are larger than the physical ones, the corresponding values 
of g^ax = (-^s — are taken as close as possible to = (see Table 1 in the next 

Section) . 

To improve the statistics, two- and three-point correlation functions have been av- 
eraged with respect to parity and charge conjugation transformations. We have chosen 
ty/a = 24 in the three-point correlation functions, which have been computed for 5 dif- 
ferent values of the initial momentum p = (27r/aL) • k, namely k = (0,0,0), (1,0,0), 
(1,1,0), (1,1,1), (2,0,0), putting always the final hadron at rest [p' = (0,0,0)]. The 
squared four-momentum transfer q'^ is thus given by = [E^^i^n) ~ ^n(E))^ ~ I^^j where 
^S(n) = ^^i(n) + kT and |gf = 0, g^n' ^Qmin^ ^iLn^ "^llin with qmin = ^Tt/uL ~ 0.7 

GeV. 

The statistical errors are evaluated using the jackknife procedure, which is adopted 
throughout this paper. 



3 Results for /i(0) 

Determination of fo{ql,ax) 

The main observation is that fo{q^) can be extracted with a statistical accuracy better 
than 0(1%) at the kinematical point = q^ax = (-^s — Afn)^ through the following 
double ratio of three-point functions with both external baryons at rest 



^ (-^ ^ N ^ Vo''^{'tx, ty, 0, 0)00 Vp^itx, ty, 0, 0)oO 

~ Vr{tx,ty,0,0)oo Vq^^ {U,ty, 0,0)oo ■ 



For large source and sink times, one has 

r> u f \ , ( (n|n7og|S-)(S-|g7on|n) \ , xi2 /^-x 

J^0[tx,ty) - — 7- — 1 \/v-|- I —[J0[qmax)\ ■ 1^0 

The double ratio (26), originally introduced in Ref. [15] for the study of heavy-heavy 
semileptonic transitions, has a number of nice features, described in detail in Ref. [1] 
and here briefly collected: 

• Normalization to unity in the SU(3) limit for every value of the lattice spacing, 
so that the deviation from one of the ratio (26) gives a direct measure of SU(3)- 
breaking effects on foiq^ax)', 
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• Large reduction of statistical uncertainties due to noise cancellation between the 
numerator and the denominator; 



• Cancellation of the dependence on the matrix elements \J Zy, and \fZ^ (see Eq. (11)) 
between the numerator and the denominator; 

• No need to improve and to renormalize the local vector current. This implies 
that discretization errors on the ratio (26) start at O(a^) and are of the form 
a^(ms — m^)^ because the ratio is symmetric under the exchange m^; 

• Quenching error is also quadratic in the SU(3)-breaking quantity (m<j — m^). 

The quality of the plateau for the double ratio (26) can be appreciated by looking at 
Fig. 1. 
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Fig. 1: Time dependence of the quantity —y/Ro{t,ty = 24a), defined in Eq. (26) for the 
two sets of the hopping parameters, given in the legend. Horizontal bars represent the 
time interval chosen for the fit. 



In Fig. 2 and Table 1 we have collected our results for the quantities /o(^maa;)> deter- 
mined through the ratio (26), the two mass combinations a'^{M^+M!^) and a'^{M^—M!^), 
the latter being proportional to the SU(3)-breaking quantity {nis — mi), and a^q'^^^ = 
a'^{MY, — Mn)"^. One can appreciate the remarkable statistical precision obtained for 
/o(9maa;)> which is always below 0.2%: the SU(3)-breaking corrections are clearly re- 
solved with respect to the statistical errors. 
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Fig. 2: Results for [— fo{Qmax)] versus the difference of the squared masses of the external 
haryons, in lattice units. 



ks - h 




a' [Ml - M^) 


^ 9max 


/o(9max) 


0.1336 - 0.1340 


0.956 (9) 


-0.0261 ( 6) 


0.000354 (17) 


-0.99995 (19) 


0.1336 - 0.1343 


0.937 (10) 


-0.0452 (12) 


0.001086 (65) 


-0.99941 (44) 


0.1336 - 0.1345 


0.924 (10) 


-0.0582 (17) 


0.00183 (12) 


-0.99858 (71) 


0.1340 - 0.1336 


0.867 (12) 


+0.0250 (8) 


0.000360 (24) 


-0.99972 (15) 


0.1340 - 0.1343 


0.824 (12) 


-0.0187 (8) 


0.000211 (20) 


-0.99985 (12) 


0.1340 - 0.1345 


0.811 (13) 


-0.0313 (14) 


0.000603 (58) 


-0.99952 (36) 


0.1343 - 0.1336 


0.783 (15) 


+0.0431 (16) 


0.001180 (99) 


-0.99933 (68) 


0.1343 - 0.1340 


0.759 (15) 


+0.0187 (8) 


0.000229 (23) 


-0.99987 (14) 


0.1343 - 0.1345 


0.728 (16) 


-0.0121 (8) 


0.000101 (14) 


-0.99992 (9) 


0.1345 - 0.1336 


0.728 (18) 


+0.0548 (27) 


0.00206 (22) 


-0.9991 (16) 


0.1345 - 0.1340 


0.704 (18) 


+0.0307 (17) 


0.000666 (84) 


-0.99977 (57) 


0.1345 - 0.1343 


0.686 (19) 


+0.0125 (8) 


0.000113 (16) 


-1.00000 (11) 



Table 1: Values of the hopping parameters ks and kc, a'^{M^ + M^), a'^(M^ — 



a {My, — Mn) and /o(9max)' obtained from the double ratio (26). 



10 



Extrapolat ion to g2 ^ 



The next step is to determine both fo{q'^) and /i(g^) for various values of q'^ using 
the standard f.f. analysis given by Eqs. (18)-(20). The quality of the plateaux for the 
quantity Wi{q'^;tx, ty) (see Eq. (16)), which provides the dominant contribution to fi{q^), 
is shown in Fig. 3. The f.f. fi{q^) turns out to be determined with quite good precision 
as it can be clearly seen in Fig. 4. Note that in Fig. 4 the data points appear paired since 
both S — > 71 and n — > S transitions are considered in the analysis. On the contrary, due 
to the large statistical noise in the quantities W2{q^] tx,ty) and Ws{q'^; tx,ty), the results 
for fo{q^) do not share the same level of precision as the one obtained for the f.f. fi{q'^), 
a finding similar to what was already observed in Ref. [1] for the K ^ it transition. 




Fig. 3: Time dependence of the quantity Wi{q'^;t,ty = 24o) (see Eq. (16)) for the two 
combinations of the hopping parameters (ksjki) = (0.1336,0.1340) (a) and {kg, kg) = 
(0.1345,0.1343) (b). The horizontal bars represent the time intervals chosen for the fit, 
which provide the value ofWi{q'^) (see Eq. (17)). Dots, squares and triangles correspond 



to = q^.^^ '^qiiin "•nd ^qi,in, with qmin = 2TT/aL, respectively. 



To improve the determination of fo{q'^) we introduce the following two double ratios 
of three-point correlation functions, corresponding to matrix elements of the spatial and 
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Fig. 4: Values of [—fi{q'^)] versus cP'c^ , as determined by Eg. (18), for the two com- 
binations of the hopping parameters {ks,ki) = (0.1336,0.1340) (a) and {ks^ki) = 
(0.1345,0.1343) (b). The solid and dashed curves represent the "dipole" and "monopole" 
fits of Eg. (30), respectively, whereas the dotted curve is a dipole fit with the slope pa- 
rameter fixed by the K* -meson mass. 



time components of the weak vector current: 

^ Im (gfclf "(tx, ty, q, 0)oo) Re (v^^{t^,ty, q, 0)oo 



R2{q'^',tx,ty) — 



Re ( Vo^''itx,ty, q, 0)oo ) Im ( qkYj^^it^, ty, q, 0)oo 



{Eg + M„)Im (gfclf ") + (Ms - MnM^Re {Vi"^) M 
{Eg + Ms)Re (V^''{t,,ty, q, 0)oo) 



Re[Vo'-^{tx,ty,q,0) 



00 



lm[qkV^^{tx,ty,q,0) 



00 



(28) 



The advantages of the ratios (28) are similar to those already pointed out for the double 
ratio (26), namely: i) a large reduction of statistical fluctuations; ii) the independence of 
the renormalization constant Zy and the improvement coefficient by ; iii) the cancellation 
of the dependence on the matrix elements \/ Zs and ^/Z^ (see Eq. (16)) between the 
numerator and the denominator, and iv) Ri{q'^]tx,ty) — > 1 in the SU(3) limit. We stress 
that the introduction of the matrix elements of degenerate transitions in Eq. (28) is 
crucial to largely reduce statistical fluctuations, because it compensates the different 
fluctuations of the matrix elements of the spatial and time components of the weak 
vector current. 
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In terms of the large-time limits Ri{q^) = lim^^ _t^) _ Ri{(f'-,tx-,ty) one has 
(Ms + Mn)h{q^) - {Eg - M„)/2(g') + {Eg + Mj:)h{q^) 



Ri{q') 



R2{q^ 



(Ms + M„)/i(g2) 
(M2 + M„)/i(g2) 



{E, 
{E, 



Mj,)h{q^) + {Eg-Mn)h{q') 
M^)f2{q^) + {Eg+Mn)h{q^) 



(Ms + M„)/i(g2) - {Eg - M^)h{q'') + {Eg - Mn)h{q') ' ^^^^ 

which can be solved in terms of /2(g^)//i(g^) and fz{q^) / fi{q^)- The latter quantities can 
be multiplied in turn by the corresponding values of fi{q^) evaluated with the standard 
procedure, obtaining in this way an improved determination of both f2{q^) and f3{q'^), 
and consequently of fo{q'^)- Our results for fo{q'^) are shown in Fig. 5, where the very 
precise lattice point at q'^ = q'^^x (^^^ rightmost one), calculated via the double ratio 
(26), is also reported. It can clearly be seen that the precision achieved for the scalar 
f.f. fo{q^) at q^ ^ Qmax is comparable to the one obtained for the vector f.f. fi{q^) for 
all the simulated quark masses. 





Fig. 5: Values of [—fo{q^)] versus aP'q^, calculated using the double ratios (28), for the 
two combinations of the hopping parameters {ksjki) = (0.1336,0.1340) (a) and {ks,k£) = 
(0.1345, 0.1343) (b). The full squares represent the very precise lattice point at q^ = q'^ax; 
calculated via the double ratio (26). The dashed and solid curves represent the monopole 
and dipole fits of Eq. (30), respectively. 

The momentum dependence of fo{q^) and fi{q'^) can be analyzed in terms of monopole 
and dipole behaviors, namely 

A 



fmon. {q ) 

fdip.{q^) 



l-g2/52 
C 

(1-^2/^2)2 



(monopole fit), 
(dipole fit). 



(30) 
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As it is clear from Figs. 4 and 5, both functional forms describe the lattice points quite 
well, with a slightly better quality in the case of the dipole form {x^ /d.o.f. ~ 0.8 1.2 
for monopole fits and /d.o.f. ~ 0.6 ^ 0.8 for dipole fits). In the case of /i(g^), 
the monopole-fit parameter B turns out to be unrelated to the value predicted by pole 
dominance (the i^*-mcson mass), while the dipole- fit parameter D agrees with the K*- 
meson mass within 15% accuracy (see dotted lines in Fig. 4). The latter finding is 
consistent with the experimental determinations of nucleon form factors. Existing data 
on the proton and neutron form factors are indeed reproduced well either by using a 
single dipole form with a radius parameter governed by the p-meson mass, or by several 
monopole terms corresponding to the dominance of the coupling of the nucleon with p, 
u> and <j) mesons and their resonances [16]. 

We have applied Eqs. (30) performing a fit to all the lattice points for both /o(g^) 
and /i(g^), and imposing the equality /o(0) = /i(0) at zero-momentum transfer. Here 
below we limit ourselves to present the results obtained in two representative cases, 
where both fo{q^) and fi{q^) are assumed to be described either by monopole forms or 
by dipole ones. Other combinations provide similar results. After fitting the parameters 
appearing in Eqs. (30) we get the values of fi{q'^ = 0) for the various combinations of 
quark masses used in the simulation. Our results are collected in Table 2 and plotted in 
Fig. 6, where they are also compared with those obtained in Ref. [1] in the case of the 
vector form factor at zero-momentum transfer, /+(0), for the K^^ decay. 



ks - ki 


a2(M|-f-M2) 


a2(M|-M2) 


/i(0) 
{monopole fit) 


/i(0) 
(dipolefit) 


0.1336-0.1340 


0.2819 (16) 


-0.01483 (12) 


-0.9970 (4) 


-0.9976 (3) 


0.1336 - 0.1343 


0.2712 (17) 


-0.02557 (23) 


-0.9905 (14) 


-0.9923 (11) 


0.133G - 0.1345 


0.2640 (17) 


-0.03271 (28) 


-0.9835 (24) 


-0.98GG (18) 


0.1340 - 0.1336 


0.2529 (17) 


-F0.01439 (17) 


-0.9966 (5) 


-0.9972 (4) 


0.1340 - 0.1343 


0.2279 (17) 


-0.01061 (15) 


-0.9979 (4) 


-0.9983 (3) 


0.1340-0.1345 


0.2209 (18) 


-0.01758 (18) 


-0.9939 (11) 


-0.9951 (9) 


0.1343 - 0.1336 


0.2209 (18) 


+0.02464 (14) 


-0.9875 (24) 


-0.9902 (18) 


0.1343 - 0.1340 


0.2068 (18) 


-h0.01050 (10) 


-0.9975 (5) 


-0.9981 (4) 


0.1343-0.1345 


0.1895 (18) 


-0.00676 (10) 


-0.9989 (3) 


-0.9991 (2) 


0.1345 -0.1336 


0.2002 (19) 


-F0.03111 (17) 


-0.9766 (60) 


-0.9818 (45) 


0.1345 - 0.1340 


0.1862 (19) 


-F0.01714 (13) 


-0.9925 (22) 


-0.9942 (17) 


0.1345 - 0.1343 


0.1759 (19) 


+0.00685 (9) 


-0.9988 (4) 


-0.9991 (3) 



Table 2: Values of the hopping parameters ks and ki, a'^{M'^ + M^), a^{Mj. - M^) and 
fi{q'^ = 0), obtained assuming either a monopole or a, dipole momentum dependence of 
both fo{q^) o.nd fi{q^) to perform the extrapolation to q^ = [see Eqs. (30) J. 

In agreement with the AdemoUo-Gatto theorem our results for hyperon semileptonic 
decays exhibit an approximate linear behavior in the quadratic SU(3)-break;ing parame- 
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Fig. 6: Results for [—fi{0)] vs. the mass difference a^{M'j^ — M^)^, obtained after ex- 
trapolation to = through either monopole (open squares) or dipole (full circles) fits 
to the lattice points for fo{q^) and fi{q^) (see Eqs. (30)). The open diamonds are the 
results obtained in Ref. [1] for the vector form factor at zero-momentum transfer, /+(0), 
corresponding to the Ki^ decay. The dotted, dashed and solid lines are linear fits to the 
data, according to the AG theorem. 



ter (M^ — M^)^ oc {nig — m^)'^, as shown in Fig. 6. Note that: i) the amount of SU(3) 
breaking in the vector form factor at zero-momentum transfer appears to be larger in 
the hyperon case with respect to the Ki^^ decay; ii) the statistical uncertainty on /i(0) 
is in the range 0.02 0.6% at the simulated quark masses, and iii) lattice artifacts 
on /i(0) due to the finiteness of the lattice spacing start at 0{a^) and are proportional 
to {nis — m^)'^, like the physical SU(3)-breaking effects; we therefore expect that, by 
having in our lattice simulation ~ 2.6 GeV, discretization errors are sensibly smaller 
than the physical SU(3)-breaking effects. Further investigations at different values of 
the lattice spacing could better clarify this point. 

Chiral extrapolation 

Following Ref. [1] we construct the AG ratio defined as 

in which the leading meson mass dependence predicted by the AG theorem is factorized 
out. This quantity depends on both kaon and pion masses, or equivalently on the two 
independent mass combinations a^(M^ + M^) and a^(M^ — M^). As shown in Fig. 7, 
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the dependence of the AG ratio (31) on the meson masses is well described by a simple 
linear fit: 

R{MK,M^) = Ro + Ri-a''{Ml + Ml) , (32) 

whereas the dependence upon the variable a^(M|^ — M^) is found to be negligible. In 
order to investigate the stability of extrapolation to the physical point we consider also 
a quadratic fit in the meson masses: 

R{Mk, M^) = Ro + Ri- a^{Ml + Ml) + R2 ■ a^{Ml + M^f . (33) 

The consequent increase of the number of parameters leads to larger uncertainties, while 
the shift in the central values remains at the level of the statistical errors. 

The extrapolation of R{Mk, M^^) to the physical point is performed using the physical 
values of meson masses in lattice units calculated by fixing the ratios Mk/Mk* and 
M^/Mp to their experimental values, obtaining 

[aMKY'^y^ = 0.189(2) , [aM^]^^^^ = 0.0536(7) . (34) 

Prom Fig. 7 it can clearly be seen that: 

• as already observed, the AG ratio for hyperon decay is ~ 2 3 times larger than 
the corresponding ratio obtained for the K^^ decay in Ref. [1]; 

• the minimum value of a^(M^ + M^) reached in our calculations is two times 
larger than the corresponding minimum value achieved in the case of the K^^ 
decay. The reason is that the quality of the signal coming from the ground state 
depends on the energy gap between the ground and the excited states (sharing the 
same quantum numbers). As the quark mass decreases, the energy gap slightly 
decreases in the case of hyperons [17], while it increases in the case of pseudoscalar 
mesons. The net result is an increase of the statistical errors at large-time distances, 
which is more pronounced in the case of hyperons with respect to the case of 
pseudoscalar mesons. Note also that we have adopted single, local interpolating 
fields for hyperons (see Eq. (6)). The use of smeared source and sink, as well as 
the use of several, independent interpolating fields may help in achieving a better 
isolation of the ground-state signal. 

Since in our simulation quark masses are rather large, we expect that the effects 
of pseudoscalar meson loops may be suppressed, which means that chiral logarithms 
are unlikely to affect significantly the lattice results. The analog case of the K n 
transition, in which the leading chiral corrections in the quenched approximation have 
been explicitly evaluated [1] , shows indeed that their contributions are almost negligible 
at the quark masses used in the simulation, as it is illustrated in Fig. 7 (compare full 
diamonds with open squares). Therefore our results for the AG ratio of the vector 
f.f. should be considered mainly as an estimate of the contributions from local terms in 
the chiral expansion. We then extrapolate the lattice results to the physical values of the 
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meson masses by assuming a simple polynomial dependence, as in Eqs. (32) and (33). 
Effects of meson loops, evaluated in full QCD at the physical quark masses, should be 
eventually added. Their impact and the actual limitations of HBChPT will be discussed 
in the next Section. We also note that, in the case of K^^ decays, preliminary results of 
unquenched simulations [2] seem to indicate the smallness of quenching effects. 



60 
40 
20 




T 1 1 1 1 1 1 1- 



(a) monopole fits in q' 



o 
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(M ^ + M ') 

K 71 

Fig. 7: Values of the AG ratio R{Mk, M^^) versus the mass combination a'^{M'^ + M^), 
for the cases of the monopole and dipole fits of both fo{q^) and fi{q^). The dashed lines 
represent the results of the linear fit (32), while the open dots indicate the value of the 
AG ratio extrapolated to the physical point. In (a) the full diamonds and the open squares 
are the corresponding results obtained for the AG ratio in the case of the K^^ decay in 
Ref. [1] with and without the subtraction of the leading chiral correction, respectively. 



In Table 3 we collect the extrapolated values R{M^^^' 
various functional forms assumed in fitting both the g^-dependence of the form factors 
(Eq. (30)) and the meson mass dependence given by Eqs. (32)-(33). Two sets of results 
are presented in Table 3. The first set, given in the upper table, is obtained by deter- 
mining the matrix elements ^/Z^ and ^/Z^ (see Eq. (16)) from a fit of the two-point 
correlation functions in the same time interval chosen for the three-point correlators, 
namely t/a £ [10,16]. For the second set, given in the lower table, the same time in- 
terval chosen to extract hadron masses, i.e. t/a £ [19,25], has been considered. Note 
that the uncertainties in the extrapolated value R{M'^^'^' , M^^^^') turn out to be quite 
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large mostly because of the long chiral extrapolation. Simulations at lower quark masses 
would be very beneficial in rochicing STich an Tinccrtainty. 





Linear Fit 


Qiiadratic Fit 


Monopole fits in 


41 ± 13 


52 ±29 


Dipole fits in 


32± 10 


35 ±30 






Linear Fit 


Quadratic Fit 


Monopole fits in (7^ 


55 ±19 


78 ±39 


Dipole fits in q^ 


40 lb 13 


49 ±35 



Table 3: Results for the AG ratio at the physical point, R{M^^^' , M^^^^') , obtained from 
linear (32) or quadratic (33) fits in the meson masses, assuming either monopole or 
dipole functional forms for the extrapolation of the scalar fo{q^) and vector fi{q^) form 
factors to q^ = 0. Upper and lower tables correspond to different choices of the time 
interval chosen for the fits of the two-point correlation functions (see text). 

From the spread of the results given in Table 3 we quote R(M^^^^' , M^^^^^' ) = 48 ± 
14stat. ± 23syst., which implies 

1 + /i (0) = RiM^""'- , MP^y'- ) • UaM^^'f - {aMP^'y'-f] ^ 

= 0.052 ± 0.015stat. ± OmSsyst. = (5.2 ± 2.9)% , (35) 

where the systematic error includes the uncertainties due to the extrapolation in q'^ and 
in the meson masses, while it does not include the uncertainty due to the quenching 
effects. In addition, the result (35) does not include the effects of chiral logarithms 
induced by meson loops, which will be discussed in the next Section. 

4 Chiral corrections in HBChPT 

The chiral behavior of the vector form factor at zero-momentum transfer /i(0) has been 
recently investigated in Ref. [7] using HBChPT, where baryons are treated as heavy 
degrees of freedom and a 1/Mb expansion around the non-relativistic limit is performed. 
The chiral corrections to /i(0) can be schematically expressed as 

.(o,./-(o) . {^-(^)-(<^1^). 

where /f^^^^(O) is the value of /i(0) in the SU(3) limit which is fixed by the vector 
current conservation. In the r.h.s. of Eq. (36) the second term in the brackets represents 
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the one-loop 0{p^) correction, while the subsequent two terms are the (parametrically) 
sub-leading 0{p^) and 0{1/Mb) corrections, respectively. The various chiral corrections 
in Eq. (36) have been computed in Ref. [7], which completes and corrects previous de- 
terminations presented in the literature [18, 19]. The corresponding numerical estimates 
are collected in Table 4 for various hyperon decays. An important feature of the chiral 
expansion (36) is that the C(p^), 0{p^) and 0{1/Mb) corrections are independent from 
unknown LECs thanks to the AG theorem. 



/i(o)//r^'^(o) 


/f"^^^(0) 


0{p^) 


0{p^) 


0{1/Mb) 


All 


^ n 


-1 


+0.7% 


+6.5% 


-3.2% 


+4.1% 


A ^ p 


-V3/2 


-9.5% 


+4.3% 


+8.0% 


+2.7% 




Y^3/2 


-6.2% 


+6.2% 


+4.3% 


+4.3% 




1/V2 


-9.2% 


+2.4% 


+7.7% 


+0.9% 



Table 4: Chiral corrections at the physical point estimated in Ref. [7] for various hyperon 
decays. 

The sum of all the contributions for /i(0)//f'^''^'*(0), presented in Table 4, turns out 
to be positive and of the order of few percent for the various hyperon decays considered. 
However the final results come from partial cancellation of larger terms. Thus higher 
order corrections are expected to give non-negligible contributions, and the convergence 
of the chiral expansion for hyperon transitions turns out to be questionable [7]. 

In addition an important source of uncertainty is represented by the mixing with 
the = 3/2+ decuplet in the effective field theory calculations. If the mass-shift A 
between the decuplet and the octet hyperons were much larger than the interaction 
scale Kqcdi the decuplet contribution could be integrated out and reabsorbed into the 
LECs, so that no correction would appear up to 0{p'^) in the chiral expansion. However 
A ~ 230 MeV is of order Aqcd and therefore the decuplet may give non-negligible, 
non-analytic contributions to the chiral expansion. 

The HBChPT with explicit decuplet degrees of freedom was firstly proposed in 
Ref. [20] and formalized as an expansion in Ref. [21]. Its impact at 0{p^), 0{p^) and 
0{1/Mb) has been investigated in Ref. [7]. As for the octet contributions, the AG theo- 
rem protects the decuplet corrections from unknown LECs and the only new parameter, 
besides A, is the known decuplet-octet-meson coupling C ~ 1.6 [22]. At 0{p^) the 
dynamical decuplet gives an important contribution to the n transition (—3.1%). 

At 0(p^) there are two contributions. The first one is due to the insertion of decuplet 
mass-shifts and it is of order —1.8%. The second contribution is due to baryon mass- 
shift insertions and gives a huge contribution of about —38%. This effect completely 
breaks the chiral expansion raising serious doubts on the consistency of the HBChPT 
with dynamical decuplet. The reason why this effect was not noticed before is because 
other quantities, at this order, contain a large number of LECs which are difficult to 
estimate from the data. In the case of /i(0) there are no LECs and a true test of the 
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convergence of the chiral expansion is possible. 

It is therefore clear that a model-independent estimate of chiral corrections for hy- 
peron decays cannot be given at present. We can limit ourselves to trust HBChPT 
without dynamical decuplet, making the Ansatz that the decuplet contributions can be 
reabsorbed into local terms. Under this assumption the chiral correction to /i(0) for 
the E~ — > n transition can be estimated from Table 4 to be (—4 ± 4)%, assuming a 
100% overall uncertainty. Thus (at least) in the case of the S~ — >^ n transition there 
seems to be a partial cancellation between chiral loop corrections from the HBChPT 
and local contributions from the lattice calculation, Eq. (35), though within quite large 
uncertainties. Adding the two results we get 

/i (0) = -0.988 ± 0.029iattice ± 0.040HBChPT . (37) 

Unquenched simulations with light quark masses are required to provide a reliable esti- 
mate of SU(3)-breaking corrections to hyperon form factors without relying on the chiral 
expansion. 

5 Results for ^/i(0)//i(0) 

The ratio of the axial to the vector f.f. at zero- momentum transfer, ^i(0)//i(0), is an 
important ingredient in the analysis of experimental data on hyperon semileptonic de- 
cays (see Refs. [3] and [23]). In this Section we first present the lattice results for the 
momentum dependence of the ratio fl'i((/^)//i(g^), then we extrapolate it in down to 

= and finally we analyze its mass dependence to obtain the ratio gi (0) / /i (0) at the 
physical quark masses. 

Using the standard procedure, the two f.f.'s gi{q^) and fi{q^) can be separately 
calculated for g / through Eqs. (23) and (18) in terms of the quantities W^j^\q^) 
and Wj{q^) = 1,2,3), defined in Eqs. (22) and (17), respectively. In this way one 
obtains values of gi{q^) with a statistical precision similar to the one achieved for the 
f.f. /i(g'). 

A reduction of the statistical noise can be obtained by considering directly the quan- 
tity gi{q^) / fi{q^) in terms of the ratios W^j^\q^) /W j{q^) . In this way the statical noise 
introduced by the two-point correlation functions in the denominators of Eqs. (21) and 
(16) is cancelled out. Moreover, also the systematic uncertainty due to different ways of 
determining the matrix elements \JZy, and v^Z^ (typically a ~ 10% effect) is removed. 
The values of gi{q^) / fi{q^) obtained at g / (i.e., q^ / l^ax) ^"^^ shown in Fig. 8 for 
two combinations of quark masses. 

At g = 0, corresponding to q^ = q^ax^ neither gi nor /i can be calculated directly. 
For both initial and final hyperons at rest there are only two non-vanishing matrix 

elements with 7 = 7' = 0, Im (^f "(t^;, ij,, (), ())oo) and Re(v^'^{t:c,ty,6,6)oo)- Their 
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Fig. 8: Values of gi{q^) / fi{q^) versus c?q^ , obtained for the two combinations of the 
hopping parameters {ks,ki) = (0.1336,0.1340) (a) and {ks,ki) = (0.1345,0.1343) (b). 
The dashed and dotted lines represent linear and quadratic fits in , respectively. 



ratio provides a new combination of f.f.'s, namely 
Im(^f-(t,,tj„0,0)oo^ - N , M. 



9\{q 



2 ) 

max) 



-92{qmax) _ giiqliax) 



Yie{V^^{t^,ty,Q,Q) 



00 



hid. 



max) I M^+Mn-' '^^^"'^"■^ ' 



k{q 



max) 
max) 



(38) 



To get giiqmax) / fiilmax) ^PP^Y ^wo Corrections to Eq. (38). The first one corresponds 
to multiply the values of giiqliax) / h{(lmax) by foiqmax)/ hiQmax), where the values 
of foiq^ax) obtained with the double ratio (26), while those for fiiq^ax) be 
evaluated using the dipole fit of Eqs. (30). The second correction is the subtraction 
of the contribution due to g2{qmax) / fii^max)i whose values can be obtained by fitting 
the momentum dependence of the ratio giiq^) / fi{q^) (see next Section). Thanks to the 
smallness of the difference My, — M„ in our simulation, both corrections turn out to be 
much smaller than the statistical errors, so that the ratio (38) can be safely corrected 
providing a determination of gi{qmax) / fiilmax) with quite good statistical precision 
(< 10%), as it is shown by the rightmost point in Fig. 8. 

The addition of a precise determination at q^ — qmax crucial to improve the 
extrapolation to q^ = 0. Indeed, thanks to the closeness of the q^ax values to q^ = 0, 
the extrapolated values g'i(0)//i(0) are only slightly affected by the specific functional 
form adopted for describing the momentum dependence of gi{q'^)/ fi{q'^), as it is can 
be clearly seen from Fig. 8. In Table 5 and in Fig. 9 we show the results obtained for 
5i(0)//i(0) using a linear fit in q"^. 

Since the axial- vector f.f. at zero- momentum transfer, 51 (0), is not protected by the 
AG theorem against first-order corrections in the SU(3)-breaking parameter (ms — m^) oc 
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kg — kf 


giVJ)/ Jiiuj 


0.1336 — 0.1340 


—0.266 (18) 


U.loDO — U.1343 


— O.zol (19) 


0.1336 — 0.1345 


—0.257 (19) 


0.1340 — 0.1336 


—0.273 (20) 


0.1340 — 0.1343 


—0.267 (22) 


0.1340 — 0.1345 


—0.264 (22) 


0.1343 — 0.1336 


—0.273 (22) 


n 1 "^zL"^ — n 1 340 


—0 975 


0.1343 -0.1345 


-0.270 (26) 


0.1345 -0.1336 


-0.274 (27) 


0.1345 - 0.1340 


-0.278 (28) 


0.1345 - 0.1343 


-0.279 (31) 



Table 5: Values of the ratio 5i(0)//i(0), obtained from a linear fit in of the lattice 
results for gi{q^) / fi{q^) ■ 

o^(M^ — M^), we start the analysis of the mass dependence of 5'i(0)//i(0) by considering 
a linear fit in the two mass combinations a^{M]^ + M^) and a^(M|- — M^): 

^ = Ao + ^1 • a\Ml + Ml) + A2 ■ a\Ml - Ml) . (39) 

The resulting values of the parameters are: Aq = —0.281 (59), Ai = 0.05 (20) and 
A2 = —0.25 (12), which provide at the physical meson masses the value 



/i(0) 



= -0.287 ± 0.052 , (40) 



where the error does not include the quenching effect. Assuming A2 = one gets Aq = 
-0.300 (61), Ai = 0.14 (21) with [gi{0)/ fi{0)]P''y'- = -0.295 ± 0.053, whereas adopting 
Ai = one gets Ao = -0.270 (21), A2 = -0.30 (26) with bi(0)//i(0)]^'''f*- = -0.280 ± 
0.026. In both cases the values extrapolated to the physical point agree with the result 
given by Eq. (40) within the statistical errors. The two separate linear fits in the mass 
combinations a^(M|. + Ml) and a^(M^ — Ml) are shown in Fig. 9 by the dashed lines. 
We have checked that the results presented for [51 (0)//i (0)]*''*^*' remain unchanged if 
instead of the meson masses the hyperon ones are adopted in the fitting procedure based 
on the Ansatz (39). We also tried quadratic fits in the mass combinations a^(M|- + M^) 
and a^(M^ — Ml). The consequent increase of the number of parameters leads to larger 
uncertainties, but the shift in the central values remains smaller than the statistical 
errors. 

The case of degenerate quark masses rUg = (i.e., q^ax = 0) allows us to access 
directly the values of 5i(0)//i(0) in the case of exact SU(3) symmetry. Without any 
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Fig. 9: Results for the ratio 5i(0)//i(0) vs. the meson mass combinations (M'^ + M^) 
(a) and a'^{M]^—M^) (b), obtained after extrapolation to = through linear fits in . 
Full dots and squares correspond to lattice results calculated with non- degenerate and de- 
generate quark masses, respectively (i.e. with and without SU (3) -breaking corrections). 
The dashed and dotted lines are linear fits in the meson masses putting in Eq. (39) 
A2 = (a) and Ai = (b), while the open dot and square are the corresponding va- 
lues extrapolated to the physical point. The diamond represents the value of the ratio 
ffi(0)//i(0) adopted in Ref. [3]. 



extrapolation in q^ we obtain the four values of (7i(0)//i(0) shown in Fig. 9(a) by the 
full dots. A linear extrapolation in the squared pion mass gives 

^^f-^ = -0.269 ± 0.047 , (41) 

LJi(0)Jst/(3) 

which nicely agrees with the prediction of the linear fit shown in Fig. 9(b) at the SU(3)- 
symmetric point a^(M^— M^) = 0. The comparison with result (40) indicates that at the 
physical point the ratio 5(i(0)//i(0) is affected by moderate SU(3)-breaking corrections, 
though it is not protected by the AG theorem. Our finding is in qualitative agreement 
with the exact SU(3)-symmetry assumption of the Cabibbo model [9]. 

Our result (40) is also consistent within the errors with the value 5i(0)//i(0) = 
—0.340 lb 0.017 adopted in the recent analysis of hyperon decays of Ref. [3], where the 
exact SU(3)-symmetry assumption for the ratio 5'i(0)//i(0) is avoided. 
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6 Results for the other form factors 



In this Section we analyze the momentum dependence of the ratios of the remaining 
f.f.'s f2{(f), 52(9^), hig'^) and gz{Q^) to the vector f.f. f\{q^). The contributions of both 
fsio^) aiid 53(9^) to the decay distributions are suppressed by the ratio of the lepton to 
the hyperon mass, and therefore they are negligible in the case of hyperon decays with 
electrons, while they can be relevant for muonic decays. The latter however suffer from 
quite small branching ratios [23]. 

As described in Section 3 the ratios f2{(f') / f and fz{<i^)/fi{(f') can be calculated 
using the double ratios (28)- (29). The latter can be easily generalized to the axial sector 
to determine the ratios g2{<f)/fi{(f) and gsiq^)/ fi{q^)- We introduce the following two 
double ratios of three-point correlation functions, corresponding to matrix elements of 
spatial and time components of the weak axial current: 



Re (^A^'^iU, ty,q, 0)00) Im (A^^{U, ty, q, 0)oo) 
Im (^A^Htx, ty, q, 0)00) Rc (^A^^iU, ty, q, 0), 



j,{A), 2. , . Re (Ag") + {Mn - Ms) Im (^f ») /q. 

Re (A^''{tx,ty,q, 0)00) 
lm(^A^^{tx,ty,q,6)oo^ 



Re(A^^it^,ty,q,6) 



00 



(42) 



We emphasize that the ratios (42) are exactly equal to unity in the SU(3) limit, and 
therefore the deviation of these ratios from one are a measure of SU(3)-breaking cor- 
rections. This is crucial for obtaining a determination of f.f.'s like the weak electricity 
g2{q^) (and the induced scalar fz{(^) in the vector case), which vanish identically in the 
SU(3) limit. 

In terms of the large-time limits R^jl^\q^) = lim^^^ _ _ ^ R\^\q'^;tx,ty) one has 
2^ (^S + M„)(/i(g2) - [Eq + Ms)<72('?') + {Eq - M^q^) 



(Ms + M„)ffi(g2) - {Eg - MMq^) + (Eg - M^q^) ' 

-p(A). 2^ ^ (Ms + M„)gi(g2) - {Eg + MMq^) + (Eg - Mj,)g^{q^) 
2 ^ (Ms + MMq") - {Eg - MMq^) + {Eg - M^q^) ' ^ ^ 

which can be can be easily solved in terms oi g2{q'^) / gi{q'^) and g's,{q^) / gi{q'^) ■ The latter 
can be multiplied in turn by the values of gi{q^) / fi{q^) to get the ratios g2{q^) / fi{q^) 
^ndg^{q^)/h{q^). 

The study of the momentum dependence of the four ratios f2{(^)/fi{(^)-, fz{q^) / fi{q^)-, 
92{q^) / fi{q'^) and g?,{q^) / f i{q^) suffers however of a crucial difference with respect to the 
case of the ratio gi{q^)/ fi{q^), namely the absence of a determination at q^ = q^^x 
in general, of a lattice point sufficiently close to = 0. This limitation makes the ex- 
trapolation to = plagued by larger uncertainties. A way to cure this problem may 
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be the use of twisted boundary conditions for the quark fields [24]. Indeed, as shown 
in Ref. [25], non-periodic boundary conditions may greatly help in determining form 
factors at zero-momentum transfer, when the latter cannot be determined sufficiently 
close to = using periodic boundary conditions. 

We discuss first the ratios f2{<l^) / fiio^) and gsiq"^)/ fi{q'^), and then the ratios 
92{<f') / fii'f') and fsiq^)/ fi{q'^), because the former have a non-vanishing SU(3) limit, 
while the latter, being due to second-class currents [26], are totally generated by SU(3)- 
breaking corrections. 

Results for /2(0)//i(0) and 53(0)//i(0) 

The momentum dependence of the ratios /2(q'^)//i(9^) and gsiq'^) / fi{q'^) is shown in 
Figs. 10 and 11, respectively. In the case of g^{q^) / fi{q^) the statistical error is very 
large for the kinematical point corresponding to initial momentum p'= {2Ti/aL) ■ (2, 0, 0) 
and therefore this point is not reported in Fig. 11. 

The absence of a determination sufficiently close to = as well as the limited num- 
ber of points introduces some sensitivity of the values extrapolated at zero-momentum 
transfer to the specific functional form assumed for the g^-dependence of the f.f.'s. The 
values of the ratios extrapolated to zero-momentum transfer adopting a linear fit in q'^ 
are collected in Table 6 and shown in Fig. 12. Note that the uncertainties are always 
larger than 50%. 



kg - ki 


/2(0)//i(0) 


53(0)//i(0) 


0.1336 - 0.1340 


-0.18 (20) 


4.9 (2.7) 


0.1336 -0.1343 


-0.20 (21) 


4.5 (2.7) 


0.1336 - 0.1345 


-0.21 (22) 


3.9 (2.6) 


0.1340 - 0.1336 


-0.31 (26) 


5.6 (2.9) 


0.1340 - 0.1343 


-0.34 (30) 


4.8 (3.1) 


0.1340 - 0.1345 


-0.36 (31) 


4.3 (2.9) 


0.1343 - 0.1336 


-0.44 (31) 


5.9 (3.2) 


0.1343 - 0.1340 


-0.46 (34) 


5.7 (3.4) 


0.1343 - 0.1345 


-0.50 (37) 


5.0 (3.5) 


0.1345 -0.1336 


-0.56 (34) 


6.0 (3.4) 


0.1345 - 0.1340 


-0.58 (37) 


5.8 (3.6) 


0.1345 - 0.1313 


-0.60 (39) 


5.5 (3.8) 



Table 6: Values of the ratios /2(0)//i(0) and g'3(0)//i(0), obtained hy performing a linear 
fit in q^. 

The mass dependence of /2(0)//i(0) is well described by a simple linear fit in a^(M^-|- 
M^), whereas the dependence upon the variable a^(M|- — M^) is found to be negligible. 
In the case of 53 (0) / /i (0) the findings are opposite and the dependence upon the variable 
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Fig. 10: Values of f2{q^) / fi{(f') versus a?q^ for the two combinations of the hopping 
parameters (A;^,^;^) = (0.1336,0.1340) ("aj and (/c^, A;^) = (0.1345, 0.1343) (h). The dashed 
and dotted lines represent linear and quadratic fits in q^ , respectively. 
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Fig. 11: The same as in Fig. 10 hut for the ratio gziq^)/ fi{q^)- 
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a2(M|, + Ml) is negligible. Using linear fits in a?{Ml. + Ml) for /2(0)//i(0) and in 
a^{Mj^-Ml) for 53(0)//i(0), we obtain at the physical point: /2(0)//i(0) = -1.14±0.66 
and 53(0)7/1 (0) = 6.3 it 3.5. In order to investigate the stability of the extrapolation we 
consider also a quadratic fit in a^{Ml^ + Ml) for /2(0)//i(0) and in o?{Ml^ - Ml) for 
(73(0)//i(0), obtaining the following results that we quote as our final estimates of these 
quantities 

-1.52 ±0.81, (44) 
+6.1 ±3.3. (45) 
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Fig. 12: Results for the ratios /2(0)//i(0) (a) and 53(0)7/1(0) (h) vs. the meson mass 
combinations a'^{Mj^ + Ml) in (a) and {M"^ — Ml) in (h), obtained after extrapolation 



to q = through linear fits in q . The dashed (dotted) lines are linear (quadratic) fits and 
the open dots (squares) are the corresponding values extrapolated to the physical point. 
The diamond represents: in (a) the experimental value of /2(0)//i(0) from Ref. [10], 
and in (b) the value 0/(73(0)7/1(0), obtained using the generalized Goldberg er-Treiman 
relation [12] and the axial Ward Identity (see text). 



In the case of the weak magnetism f.f. the experimental result for /2(0)7/i(0) is 
known from Ref. [10], namely /2(0)7/i(0) = —1.71 ± 0.12stat. i 0.23syst.5 which is shown 
in Fig. 12(a). In the case of 53(0)7/1(0) instead there is no direct experimental informa- 
tion. By combining the generalized Goldberger-Trieman relation [12] and the axial Ward 
identity one can argue that the f.f. 53 (g^) should have a if-meson pole at low values of 
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g^, which imphes that close to the chiral hmit one has 

93(0) ^ f M^ + Mn \ ' 5i(0) 

/i(0) \ Mk J /i(0) • ^ ' 

Using the result (40) we get 53(0)//i(0) = 5.5 it 0.9, which is consistent with Eq. (44). 
This estimate is also presented in Fig. 12(b). Note that both the latter and the experi- 
mental value for /2 (0) / /i (0) are quite closer to the results (45) and (44) , respectively. 

A calculation of the weak magnetism /2 and the induced pseudoscalar gs f.f.'s for 
values of closer to = and at lower values of the quark masses is mandatory to 
improve the precision of the lattice determination of these f.f.'s. 



Results for 92(0) /fi{0) and /3(0)//i(0) 

The momentum dependence of the ratios g2{Q^) / fi{<f) and /3((?^)//i((?^) is presented 
in Figs. 13 and 14, respectively. In the SU(3) limit both the weak electricity §2(0.^) and 
the induced scalar fz{q^) f.f.'s vanish identically, so that the results shown in Figs. 13 
and 14 arc totally generated by the breaking of the SU(3) symmetry. As in the case 
of g^iff) / fii'f') the statistical error for g2{(f')/ fi{(f') becomes very large for the largest 
absolute value, which therefore is not presented in Fig. 13. Note that in Figs. 13 and 
14 we have explicitly reversed the signs of the values of 52(9^) and fz{(]^) corresponding 
to the n ^ E transitions, because both (72(9^) and f-^{q^) belong to second-class currents. 

The sensitivity of the extrapolated values at zero-momentum transfer to the spe- 
cific functional form assumed for the g^-dependence of the two ratios 5'2(<?^)//i(<?^) and 
fz{(^)/fi{q^) is much more limited with respect to what has been observed in the case 
of the ratios f2{q^)/h{q^) and g^{q^) / h{q^) . The values of ff2(0)//i(0) and /3(0)//i(0), 
obtained adopting a linear fit in (7^, are presented in Table 7 and in Figs. 15(a) and 16(a). 
Note again that, as expected for second-class currents, the signs of both 5'2(0)//i(0) and 
/3(0)//i(0) are linked to the sign of a^iM]^ -Ml). 

Since the values of 52(0)7/1(0) and /3(0)//i(0) are exactly known in the SU(3) limit 
(i.e., 5'2(0)//i(0) = /3(0)//i(0) = 0), the analysis of the mass dependence of these ratios 
can proceed in a way similar to the one adopted for /i(0). In this case, however, taking 
into account that both g2 and are not protected by the AG theorem, we introduce 
the following ratios: 

The mass dependence of these quantities is well described by simple linear fits in the mass 
combination a^(M^ -|- M^), whereas the dependence upon the other variable a^(M^ — 
M^) turns out to be negligible. At the physical point we get Rg^{M^^'^' ,M^^^^''') = 
18.0±7.5andi?/3(M^''^*-,M^''^* ) = -11.9±6.2. A quadratic fit in the mass combination 
a'^{M'^ + Ml) does not modify the central value of Rg^, but it increases the uncertainty, 
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Fig. 13: Values of g2{q^) / fi{(f') versus a^q^ for the two combinations of the hopping 
parameters (A;^,^;^) = (0.1336,0.1340) ("aj and (/c^, A;^) = (0.1345, 0.1343) (h). The dashed 
and dotted lines represent linear and quadratic fits in q^ , respectively. 
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Fig. 14: The same as in Fig. 13 but for the ratio fz{q^) / fi{q^ 
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(J. 133d — U.134U 


—0.083 (38) 


+0.065 (14) 


L).1o3d — U.1343 


—0.139 (70j 


+0.123 (28) 


0.1336 — 0.1345 


—0.185 (92) 


+0.166 (40) 


0.134U — 0.1335 


+0.100 (43) 


r\ r\'i'\ ( 1 (W 

— 0.071 (19) 


0.134U — 0.1343 


n no r / o 'i\ 

—0.085 (37) 


+0.064 (19) 


0.1340 — 0.1345 


—0.132 (66) 


+0.112 (35) 


0.1343 — 0.133o 


+0.2z0 (88) 


C\ 'X Af* (A T\ 

—0.146 (47) 


0.1343 — 0.1340 


1 f\ C\C\ A (A C\\ 

+0.094 (42) 


—0.070 (25) 


n 1 QAQ _ n 1 S45 


— n 071 

VJ.VJ 1 ± I OZ; / 




0.1345 - 0.1336 


+0.326 (123) 


-0.213 (80) 


0.1345 - 0.1340 


+0.180 (77) 


-0.132 (56) 


0.1345 - 0.1343 


+0.066 (33) 


-0.057 (27) 



Table 7: Values of the ratios 52(0)//i(0) and /3(0)//i(0), obtained by performing a linear 
fit in . 



i.e. Rg^iMp'- ,Mfy'- 

larger uncertainty is found in the case of Rj^, namely Rf.^{M'^ 
Thus, using the values of Rg^ and Rf^ from the linear fits, we derive our final estimates 
of 52(0)//i(0) and fz{^)l f\{G) extrapolated to the physical point 



18.3 ± 14.6. A slight shift in the central value within a much 

-18.4+15.5. 



phys. j^phys 



52(0) 1"'^" 
/i(0). 

f^^vhys. 
L/i(0). 



+0.63 ± 0.26 , 



= -0.42 ± 0.22 . 



(49) 
(50) 



Our results indicate a positive, non- vanishing value of <j'2(0)//i(0) due to SU(3)- 
breaking corrections. In the experiments neither S'i(0)//i(0) nor 52(0)//! (0) are sepa- 
rately determined, but only a specific combination of these ratios can be extracted. The 
following combination has been determined for the S~ — > n transition [10] 



5i(0)- 0.133 -52(0) 



/i(0) 



exp. 



= 0.327 ± 0.007stat. ± 0.019syst. 



(51) 



Using our results (40) and (49) we get 

5i(0)- 0.133- 52(0) IP'^^^ 



/i(0) 



0.37 + 0.08 



(52) 



in good agreement with the experimental value. Even though the experimental data 
are compatible with general fitting procedures which make the conventional assumption 
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Fig. 15: Results for the ratios 5'2(0)//i(0) (a) and Rg^ (h) vs. the meson mass combi- 
nations a'^{M^ — M^) in (a) and a'^{M'^ + M^) in (b), obtained after extrapolation to 
= through linear fits in . The dashed line is the result of a linear fit and the open 
dot is the corresponding value extrapolated to the physical point. 
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Fig. 16: The same as in Fig. 15 but for the ratios /3(0)//i(0) (a) and Rj^ (b). 
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52 ('7^) = (as done for instance in Ref. [3]), it has been found that positive values for 
52(0)//i(0) combined with a corresponding reduced value for |g'i(0)//i(0)| are slightly 
preferred [10]. The lattice results (40) and (49) appear to favor the second scenario. 

7 Conclusions 

We have presented a lattice QCD study of SU(3)-breaking corrections in the vector and 
axial form factors relevant for the hyperon semileptonic decay S~ n I v. Though our 
simulation has been carried out in the quenched approximation, our results represent the 
first attempt to evaluate hyperon form factors using a non-perturbative method based 
only on QCD. 

For each form factor we have studied its momentum and mass dependencies, obtain- 
ing its value extrapolated at zero-momentum transfer and at the physical point. Our 
final results are collected in Table 8, where the errors do not include the quenching effect. 



/i(0) 


-0.988 ± 0.029iattice ± 0.040HBChPT 


,9i(0)//i(0) 


-0.287 ±0.052 


,/2(0)/,/l((J) 


-1.52 ± 0.81 


/3(0)//l(0) 


-0.42 ±0.22 


52(0)//i(0) 


+0.63 ± 0.26 


53(0)//l(0) 


+6.1 ±3.3 



Table 8: Results of our lattice calculations of the vector and axial form factors for the 
S~" — >■ n transition. 

We conclude with few additional comments: 

• The SU(3)-breaking corrections to the vector form factor at zero-momentum trans- 
fer, /i(0), have been determined with great statistical accuracy in the regime of 
the simulated quark masses, which correspond to pion masses above 0.7 GeV. The 
magnitude of the errors reported in Table 8 is mainly due to the chiral extrapola- 
tion and to the poor convergence of the Heavy Baryon Chiral Perturbation Theory 
[7]. Though within large errors the central value of /i(0) arises from a partial 
cancellation between the contributions of local terms, evaluated on the lattice, 
and chiral loops. This may indicate that SU(3)-breaking corrections on /i(0) are 
moderate (at least for the transition considered) , giving support to the analysis of 
Ref. [3]. 

• The ratio gi (0) / /i (0) is found to be negative and consistent with the value adopted 
in the analysis of Ref. [3] . The study of the degenerate transitions also allowed us 
to determine the value of (/i(0)//i(0) directly in the limit of exact SU(3) symmetry, 
obtaining [gi{^) / fi{^)]su{d,) = -0.269 + 0.047. This means that SU(3)-breaking 
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corrections are moderate also on this ratio, though the latter is not protected by 
the Ademollo-Gatto theorem against fist-order corrections. 

• The weak electricity form factor at zero- momentum transfer 52(0) is found to be 
non- vanishing because of SU(3)-breaking corrections. Our result for gi{0) / fi{Q) 
combined with that of 52(0)7/1(0) are nicely consistent with the experimental 
result from [10]. Our findings favor the scenario in which 32(0)//i(0) is large and 

positive with a corresponding reduced value for |(7i(0)//i(0)| with respect to the 
conventional assumption 52 (<?^) = (done for instance in Ref. [3]) based on exact 
SU(3) symmetry. 

Finally, wc discuss few possible improvements for future lattice QCD studies of the 
hyperon semileptonic transitions: 

• The quenched approximation should be removed and the simulated quark masses 
should be lowered as much as possible in order to reduce the impact of the chiral 
extrapolation. 

• The accuracy of the ratios /2(0)//i(0), /3(0)//i(0), 52(0)//i(0) and 93(0) /fi{0) 
can be improved by implementing twisted boundary conditions for the quark fields. 
In this way values of the momentum transfer closer to = can be accessed in 
the simulations. 

• The use of smeared source and sink for the interpolating fields as well as the use 

of several, independent interpolating fields may help in increasing the overlap with 
the ground-state signal, particularly at low values of the quark masses. 
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